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Abstract 

We show that the threshold of complete synchronization in a lattice of cou- 
pled non-smooth chaotic maps is determined by linear stability along the 
directions transversal to the synchronization subspace. We examine care- 
fully the sychronization time and show that a inadequate observation of the 
system evolution leads to wrong results. We present both careful numerical 
experiments and a rigorous mathematical explanation confirming this fact, 
allowing for a generalization involving hyperbolic coupled map lattices. 
Keywords: coupled map lattices, long-range interactions, synchronization 
time 
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The possibility of synchronizing chaotic dynamics has been harnessed in 
a large number of systems of physical interest l2| , like coupled Josephson 
junctions [3] and lasers Although there have been identified different 
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4 types of chaos synchronization, we shall concentrate on the so-called ampli- 

5 tude or complete synchronization, for which all dynamical variables undergo 

6 the same time evolution [5]. The essential dynamics involved in the process 

7 of chaos synchronization lies on the low- dimensionality of the subspace (in 

8 the phase space of the system) in which synchronized motion sets in. 

9 For example, if we consider a lattice of coupled oscillators, each of 

10 them represented by a vector field of D dimensions, where typically D N, 

11 the synchronized state belongs to a D-dimensional subspace of the ND- 

12 dimensional phase space. In order for this synchronized state to exist the 

13 coupling among oscillators takes on a suitable form 6|. Whether or not this 

14 synchronized state is stable, however, is a more difficult question, since it in- 

15 volves the analysis of infinitesimal displacements from the synchronized state 

16 along all (N ~1)D directions transversal to the synchronization subspace 

17 The stability condition of the synchronized orbit with respect to transversal 
IB perturbations can be obtained from the negativeness of the largest transversal 

19 Lyapunov exponent. 

20 In this paper we consider a coupled chaotic map lattice (CML) in which 

21 the coupling prescription is linear and non-local, for it takes into account the 

22 distance between maps along the lattice. Such non-local couplings appear in 

23 many problems of physical 8| and biological interest jo]. We suppose that the 

24 coupling strength decreases with the lattice distance as a power-law, which 

25 characteristic exponent can take on any non-negative value jlo|. The loss 

26 of transversal stability of the synchronized state, as the coupling parameters 

27 are varied, was found in such power-law couplings, with help of the largest 

28 transversal Lyapunov exponent, for a number of chaotic maps ll|, |l2|. In 
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29 the particular case of maps with constant eigenvalues of the Jacobian matrix 

30 (piecewise-linear chaotic maps) we obtained analytical results for the loss of 

31 transversal stability of the synchronized state which agree with the numer- 

32 ical simulations [13]. Such CML's represent hyperbolic dynamical systems 

33 (see text below), what enables us to use powerful mathematical tools like 



34 ergodicity and global shadowing of numerically generated orbits [14 1. 

35 On the other hand, in a recent paper there was argued that in the special 

36 case of coupled non-smooth discontinuous maps the synchronization transi- 

37 tion would not be given by the largest transversal exponent, but rather by 

38 a different approach taking into account finite distances from the synchro- 



39 nized state 15|. To investigate this apparent contradiction we considered in 

40 this paper the transient behavior of the non-synchronized orbits for coupled 

41 piecewise linear maps. Our results show that the analytical results of Ref. 

42 I3I (using linear transversal stability of the synchronized state) hold for both 

43 smooth and non-smooth maps, the numerical results being strongly affected 

44 by many factors as the large transient time and the choice of initial condi- 

45 tions. Due to these factors, the time it takes to achieve convergence to the 

46 synchronized state may be extremely large, what may lead to wrong conclu- 

47 sions about the stationary state of the system. Motivated by this problem, 

48 we investigated the validity of the transversal linear stability analysis in a 

49 class of hyperbolic CML's, using periodic-orbit theory to unveil the role of 



50 the unstable orbits embedded in the synchronized state |16l4l8|. 

51 The CML we consider in this work can be written in the explicit form of 

52 a A^- dimensional dynamical system 

x„+i = (1 + C)F(x„) = BF(x„), (1) 
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53 where the components of x„ = {xn \ Xn \ ■ ■ ■ , xl^"^) denote the state vari- 

54 able attached to the map located at the site i = 1, . . . ,N at time n = 

55 0, 1, . . .. If the uncoupled maps are written as a; i— ?■ f{x) we can write 

56 F(x) = [f{x^^^),f{x^'^^),...,f{x^'^^)]^. Moreover, the coupling prescription 

57 is represented by the matrix C, and 1 is the identity matrix. 

58 In the following we consider the generalized Bernoulli map f{x) = (3x, 

59 {mod 1), where x E [0, 1) and (3 > 1, such that the isolated map generates 

60 a strongly chaotic orbit. When these piecewise-linear maps are coupled ac- 

61 cording to Eq. ([1]), in order to ensure that Xn^ G [0, 1), the elements of B 

62 must satisfy the following necessary and sufficient conditions: Bij > 0, and 

63 < XljLi < 1, for all i,j = 1,2, ... N 13]. Moreover, we use a symmetric 

64 coupling matrix with elements 

a,=£7?-^[r,7(l-5,,)-r/5,,], (2) 

65 where rij = mmi(zz\i — j + IN\ is the minimum lattice distance between 

66 the sites i and j (with periodic boundary conditions), r] = r~", with 

67 A^' = (A^ — l)/2 (which requires odd), and the coupling strength satisfies 

68 < £ < 1 due to the constraints on B^j. The effective range of interactions 

69 is represented by a > such that the limits a = and a — t- oo correspond, 

70 respectively, to global (mean field) and local (first neighbors) coupling pre- 

71 scriptions. 

A completely synchronized state is the chaotic orbit for which xil^ = . . . = 
x^n\ and which is a solution of Eq. ([1]). Since the Jacobian DF = /3B is a 
circulant matrix, its eigenvalues can be analytically obtained as 

AW = /3[(l-£) + (£/,^)6('=-i)], 
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where 



TV' 



m=l 



/ 271 km\ 

[it) 



{0<k<N) 



(3) 



such that the Lyapunov spectrum {Aj}^^ can be derived {l3|. The stabihty 
threshold of the synchronized state with respect to infinitesimal transver- 
sal displacements, obtained by imposing A2 = 0, gives two curves in the 
parameter plane {e versus a): (i) e'^{a,N) = mm{e up{a, N),!}; and (ii) 
ec(«, N) = min{e;o(Q;, A^), 1}, where we defined 



'Up 



a,N) = + 



(4) 
(5) 



In order to check the validity of these analytical conditions for the thresh- 
old of transversal stability we have made careful numerical experiments using 
the same criteria as proposed in Ref. 15| (where it has been claimed that 
those conditions would hold only for coupled continuous maps). According 
we choose initial conditions x^^ uniformly distributed in the interval [0, 1) 
The CML is firstly iterated for a transient time of = IC" x times and 
further iterated by more T = 10^ x A^ times. As a numerical diagnostic of 
complete synchronization we computed the following quantity 



R 



E- 



NT 



X. 



(6) 



which is essentially a mean deviation from the lattice-averaged amplitude. 
The resulting dynamical state is considered as being completely synchro- 
nized if -R < 10^^. In the coupling parameter space we keep a constant 



^We initialized the random number generator rani from Ref. [20[ with always the same 
the seed (-28937104) 
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89 and sweep through the values of e G [0, 1]. The value corresponding to the 

90 synchronization threshold, denoted as Enum, is obtained from bisection as 

91 £num = i^s + £^d)/2, where Eg and Ed are, respectively, the last value corre- 

92 sponding to a synchronized state and the first value for a non-synchronized 

93 one. The numerical value of Enum is turned more accurate from refining the 

94 increment mesh and repeating the process, until {Es — Ed) < 10^^. 

95 The results of this numerical procedure, for the case (3 = 1.1, are de- 

96 picted in Figure [H where we show the value of the coupling strength at the 

97 synchronization threshold as a function of a. In Fig. [I](a) we show how the 

98 numerically determined critical value increases with a for different lattice 

99 sizes A^, the transient time being different for each choice, using w = 5. The 

100 solid lines correspond to the analytical condition derived in Ref. 13j (and that 

101 depend on the lattice size as well). In fact, as the lattice size increases, 

102 the numerical values of Enum may no longer match the analytically predicted 

103 values, if a is large enough ^^]. This does not mean, however, that the an- 

104 alytical value of Ec is not valid in those cases, but rather that the numerical 
106 simulations have not been performed using a transient long enough. To show 

106 the influence of the transient time in the results, we show in Fig. [^b) the 

107 dependence of Enum with a for a fixed lattice of = 129 sites by changing 

108 the parameter w. By increasing the transient time the numerically obtained 

109 values for the synchronization threshold agree better with those derived from 
no transversal linear stability. The same conclusions were obtained using other 

111 lattice sizes as well. These results suggest that the analytical result for Ec 

112 remains valid, as long as we use sufficiently long transient times, in contrast 

113 with Ref. \m. 
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114 Another factor that affects the accuracy of numerical results for the 

115 threshold of synchronization is that a distribution of initial conditions over 

116 the interval [0, 1) should respect the natural measure of the chaotic orbit, 

117 because we are assuming the coupling between typical oscillators, which are 

118 characterized by trajectories in the steady-state system, i. e., trajectories 

119 that satisfy the invariant density of the system. While for integer values of 

120 (3 the natural measure is uniform, this is no longer valid for fractional (3, 

121 and small errors may be introduced if we choose initial conditions with a 

122 uniform probability distribution. In order to overcome this problem we iter- 

123 ated each map s times before starting coupling them according to Eq. ([1]) 

124 (this transient time should not be confused with the transient time we 

125 compute after having started coupling the maps). In Figure [T|^c) we com- 

126 pare the results of two simulations: for the line with filled triangles we used 

127 initial conditions uniformly distributed along [0,1), without discarding any 

128 transients (s = 0); whereas the line with open triangles was obtained from 

129 initial conditions chosen with respect to a numerical approximation of the 

130 natural measure, the latter having being obtained from a transient time of 

131 s = 10^ iterations. The results obeying the natural measure of the uncoupled 

132 oscillators are more likely to agree with the analytical results since, after the 

133 synchronized state sets in, the corresponding orbit must follow this natural 

134 measure. As will be formally discussed below, in the limit n — )■ oo, the results 

135 are independet of the (typical) initial distribution of trajectories. But for fi- 

136 nite time intervals, time synchronization can depend on such distribution. 

137 This dependence is due to local dynamics and form of coupling. 

138 In order to analyze the dependence of Enum on the initial conditions, we 
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139 performed extensive numerical simulations with an ensemble of 5000 identi- 

140 cally prepared CML's, each of them with a different initial condition [Fig. 

141 [It^d)] and the same transient time {w = 5). We observed different values of 
i« ^num for each initial condition, provided a is large enough. Instead of showing 

143 each of them (what would turn the figure too much loaded with symbols) we 

144 represented in Fig. [T]^d) only those numerical values of Enum that are closest 

145 (open circles) and farthest (filled circles) with respect to the analytical value 

146 (full line). Note that synchronization of one typical trajectory implies global 

147 stability of the synchronized state, since the system is hyperbolic. A chaotic 

148 invariant set fl is hyperbolic if the following conditions are fulfilled: (i) the 

149 tangent space at each point x G i7 can be decomposed in two invariant sub- 

150 spaces (a stable and an unstable one) with constant dimensions; (ii) these 

151 subspaces always intersect transversely (i.e., they cannot present tangencies); 

152 and (in) this decomposition is consistent under the dynamics in Q generated 



by F 



14( 1 . For coupled generalized Bernoulli maps the set Q is the A^-torus 



154 [0, 1)^ and the Jacobian matrix DF has constant entries and does not de- 

155 pend on X G thus the dimension of the invariant subspaces is constant 

156 everywhere [condition (i)]. Thanks to this particular form of the Jacobian its 

157 eigenvectors (which span the invariant subspaces) are everywhere orthogonal 

158 [condition (ii)]. Let u be any of such eigenvectors: under the dynamics of F 

159 it follows that u is mapped to a vector along the same direction [condition 

160 (iii)]. Hence the set is a hyperbolic structure for F. 

It is possible to understand, from a general point of view, the causes of the 
strong dependence of the synchronization threshold results on the transient 
time and the initial conditions. These causes are not restricted to coupled 
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piecewise-linear maps as ours, but are rather generic for hyperbolic CML's. 
We can extend our conclusions to a CML given by Eq. ([1]) where the coupling 
prescription keeps invariant the phase space Q = [0, 1)^, and for which 

5 = {xG^]:a;« = ■■■ = x(^)} 

161 is the one-dimensional invariant synchronization manifold defined by the cor- 

162 responding state. We consider a A-neighborhood of S as the set of points 

163 whose distances from the S do not exceed A: Sa = {x : d(x, 5) < A}, 

164 where d is a suitably defined distance on the metric space Q. We define 

165 S = limA->o as a linear neighborhood of S. Accordingly F = — S is the 

166 phase space region, except the linear neighborhood of the synchronization 

167 manifold. 

168 We can speak of the global dynamics generated by the coupled map lattice 

169 x„+i = BF(x„) in terms of their periodic points. In this spirit we denote 

170 Xj(j9) the jth fixed point of the p-times iterated vector function BF(x„). The 

171 ith eigenvalue of the Jacobian matrix of BFf^^(x„), evaluated at this point, 

172 is written as Aj(xj(p)), such that |Ai(xj(p))| > ■ ■ • > |AAr(xj(p))|. 

173 Let us consider a subset of the phase space, A G Q, with natural measure 

174 fi{A). Note that, by construction, we have = 1. For hyperbolic systems 

the natural measure of such subset can be obtained 



176 from the unstable periodic points embedded in it as 



3 



/i(A)= lim5^1/L,(p), (7) 



hyperbolic system satisfying Axiom-A must be also mixing. This condition is ful- 
filled if the system possesses a dense set of unstable periodic orbits embedded in the phase 
space [l^. 
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177 where Lj{p) = Ylili ('^m is the largest integer such that \Ad^{-Xj{p))\ > 

178 1) and the sum sweeps over all Xj(p) G A. The exploitation of this iden- 

179 tity is the object of periodic-orbit theory, that has been used for a number 

180 of theoretical investigations on the properties of chaotic dynamical systems 

181 [isj]. For generalized Bernoulli maps /3x (mod 1) and a linear coupling, 

182 the Jacobian matrix has constant entries and thus do not depend on the 

183 orbit points, i.e., all the unstable periodic orbits have the same eigenvalue 

184 spectra (consequently Lj{p) = L{p) for all j), and the natural measure is 

185 yu(A) = limp^oQ Na{p)/ L{p), where Na{p) is the number of period-p points 

186 contained in the subset A of Q. 

187 A byproduct of the periodic-orbit theory is that the (linear) transversal 

188 stability of the synchronization manifold can be studied either from the nat- 

189 ural measure of a typical chaotic orbit (by the second largest Lyapunov expo- 

190 nent) or from the atypical measure generated by the unstable periodic orbits. 

191 In particular, with respect to the period-p orbit the threshold of transversal 

192 stability of the synchronization manifold can be obtained from the condition 

193 |A2(xj(p))| = 1 for all Xj(p) G S. As the period p goes to infinity we expect 

194 an increasingly better agreement of this result with that obtained by using 
196 the second largest transversal Lyapunov exponent (or A2 = 0). For a given 

196 a and values of the coupling strengths such that eio{a) < e < Eupia), the 

197 natural measure of the subset A is 

/i(A) = hm NAip)/l3P. (8) 

p— >oo 

198 Taking A to be the linear neighborhood of the synchronization manifold, 

199 S, there follows that the number of orbits in this neighborhood is N^, = P'^ — l 

200 for integer /3 (if /3 is fractional, as in the numerical simulations of the previous 
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201 section, Np — > for p ^ 1) and the corresponding natural measure is given 

202 by 

/.(S) = lim - l)/3-P = 1, (9) 

203 demonstrating that the linear neighborhood of the synchronization manifold 

204 S is the asymptotic state of any typical initial condition (in the sense that the 

205 set of initial conditions that do not converge to E has zero Lcbesgue measure). 

206 This result is obtained for the parameter regime in which the synchronized 

207 state is locally stable and, therefore, any trajectory in E converges exponen- 

208 tially to 5 at a rate A2 < 0. An immediate consequence of this result is that 

209 the natural measure outside the linear neighborhood is zero since, using the 

210 fact that the natural measure is ergodic, wc have /x(r) = //(fi) — Ai(S) = 0. 

211 Given that almost all initial conditions outside the synchronization man- 

212 ifold eventually asymptote to it, we may well ask why sometimes it takes so 

213 long for this convergence to be observed in numerical experiments. As we 

214 saw previously, this long transient time may even be mistaken as a effect 

215 of non-convergence. The answer lies in the properties of the horseshoe-like 

216 invariant chaotic set embedded in F. This set is non-attracting since almost 

217 all initial conditions in F converge to S as the time goes to infinity. 

218 Let Pn(x) be the density of trajectories arounf x at time n, so that 

li{A) = lim / p„(x)dx (10) 

219 for any typical po(x). The above result is independent of the specific form 

220 of po in the limit n — > 00. The possibility of expanding po in terms of the 

221 eigenfunctions of the Perron-Probenius operator justifies such independence 

222 of /i(A), since the invariant density p(x) of the system is associated with the 
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223 largest eigenvalue (which is not degenerate) of that operator [21|. Note that 

224 for finite time intervals, the convergence of the invariant density depends on 

225 the coefficients of expansion of po(x) on the basis of eigenfunctions of the 

226 Perron- Frobenius operator. The results in Figure Wic) are evidence of the 

227 assertion of the previous sentence. 

228 However, the measure generated by chaotic orbits whose initial conditions 

229 are uniformly distributed over of an open region B of the phase space Q 

230 decays exponentially with time with escape rate 7, 

/ p„+^(x)dx = e"'^" / p„(x)dx (11) 
Jb J b 

231 with m ^ 1. For hyperbolic systems it can be shown that the escape rate is 

232 also obtained in terms of unstable periodic points of the saddle according to 

hme-^^ 1/^.(P) = 1, (12) 

Xj(p)GB 

233 where in the sum we consider only the periodic orbits of the horseshoe-like 

n 

234 set B outside the synchronization manifold 16[. Hence, if one picks up at 

235 random an initial condition off the synchronization manifold, the distribu- 

236 tion of the transient times is likely to be exponential, with a characteristic 

237 exponent dependent on the escape rate 7. This is illustrated in Figure [21 

238 in which the distribution of synchronization times (transient time intervals), 

239 (^in)^ is indicated for a typical realization of the network, with N = 17. 

240 Figure [2] also points out a numerical estimate of the density measure de- 

241 cay of the unsynchronized state. To obtain this estimate, we cover F with 

242 Kq uniformly distributed initial conditions and, at each time instant n, we 

243 count the number Kn of trajectories that remain in F under the system 
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244 evolution. The temporal decay of Kn provides the escape rate of T since 

245 Kn/Ko ^ /pp„(x)dx, for Kq ^ 1. Note that the exponential decay of the 

246 curves is the same in Figure [2l The mean time of synchronization (or the 

247 mean life of the chaotic saddle), given by 

{n) = / n(p{n)dn ^ — , (13) 
Jo 7 

248 is indicated by a red arrow on the x-axis. The rightmost term in Eq. ( !T3|) 

249 is obtained by supposing 0(n) oc e~^^ - which is typical for chaotic saddles 

250 and is verified in Fig. [2J 

251 If the initial condition is too close to an unstable periodic orbit (or its 

252 stable manifold) of B it would stick to it for some time-span and hence it 

253 takes a very long time for such a trajectory to approach the synchronization 

254 manifold, as illustrated in Fig. [31 This seems to occur very often if we use 

255 fractional values of /3, like in the numerical simulations we shown in this 

256 paper. 

257 We show in Figure [3] a situation in which the trajectory escapes from 

258 the neighborhood of an orbit with small period (~ 30), and instantly ac- 

259 cess the synchronized state. However, since there is a chaotic saddle in F, 

260 there is a dense set of unstable periodic points in the unsynchronized state, 

261 and a trajectory can wander between different UPO's before escaping to the 

262 synchronized state. 

263 To demonstrate that the behavior shown in Figure [2] is typical, we present 

264 in Figure H] two numerical estimates for the invariant density of the chaotic 

265 saddle in the unsynchronized state. In Figure |l](a), which is obtained from 

266 typical trajectories, the invariant density is shown in a projection on the 

267 x*^*) X plane (due to symmetry the network, the value of i is irrelevant). 
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268 The most visited regions in that figure are represented by symbols whose 

269 color is dark blue and the least visited by light green symbols. The estimate 

270 for the density in terms of UPO's is shown in Figure Hl^b). For this case, all 

271 the unstable orbits with the period less than or equal to 35 are considered. 

272 Comparing the two figures is apparent that the structures in Figure ID^a) 

273 are supported by the orbits in Figure 111(b). The processing time was the 

274 limiting factor for choosing the value of the period in the simulations shown 

275 in Fig. m^bjfl. The Table [T] shows the mean differences between the values 

276 of projections of p(x) and Pp(x), the density values obtained by calculation 

277 of all points in F with period q < p. This calculation was done in a grid of 

278 128 X 128 boxes, each one as the same size. We also found that the average 



279 difference decreases as the period increases, this is according to reference 



22| 



280 When the period tends to infinity, we obtain the exact measure of the saddle 



282 From the view of the structure of the chaotic saddle provided by Figure HJ 

283 it is easy to understand how the distribution of initial conditions affects the 

284 results for simulations with finite time. For the case we are considering, the 

285 Bernoulli map with /3 = 1.10, the density of the natural measure in the unit 

286 interval is mainly concentrated near zero. In Figure IH^b) the chaotic saddle 

287 is less dense at the bottom left than in the upper right corner, i.e., the initial 

288 conditions close to zero spend less time to synchronize. However, it should 

289 be noted that this dependence with the initial distribution of trajectories is a 

290 consequence of the finite time simulations because, as shown, such as effects 



■^For p = 28 there are around 10000 periodic points in F. However, for p = 35 there are 
around 10000000 points. 
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Period p 


(p(x)-p,(x))xl04 


29 


1.164 


30 


1.157 


31 


1.153 


32 


0.690 


33 


0.499 


34 


0.458 


35 


0.425 



Table 1: The mean difference (/9(x) — pp(x)), magnified by a 10'* factor, in a grid of 
128 X 128 boxes, each one as the same size. 



are transient due to the existence of the chaotic saddle in F. Figure H] shows 
a high value for the invariant density for specific points along the diagonal 
that contains the projection of the synchronized state. These points, which 
are in the chaotic saddle, are those associated with the discontinuity of the 
local map, and the respective pre-images. Thus, the probability of a 
trajectory in the neighborhood of such points still belongs to the chaotic 
saddle is very high. Consequently, the behavior illustrated in Figure Et^a) is 
more likely to occur when a trajectory close to S is near of a discontinuity (or 
its pre-images) of the local dynamics. This is the topological explanation for 
the non-local instabilities in synchronized state, which supports the heuristic 
argument presented in Reference [l5|. Direct analysis of our results, espe- 
cially those presented in Figures |3] and HI shows that the divergence between 
a typical trajectory and locally stable synchronized state is due strictly to 
the existence of UPO's in the unsynchronized state. 
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305 In conclusion, the analytical conditions for the threshold of transversal 

306 stability of the synchronized state of coupled piecewise-linear maps are con- 

307 firmed by numerical experiments as long as we observe the following precau- 

308 tions: (i) the transient time should be chosen as large as possible, (ii) the 

309 choice of initial conditions should be done using a probability distribution 

310 which best matches the natural measure of the uncoupled oscillators. 

311 Although these computational problems are less likely to occur in cou- 

312 pled smooth maps, they do not invalidate the analytical approach to the 

313 transversal stability of coupled non-smooth maps, like piecewise-linear ones. 

314 Our analysis indicates that results based on linear analysis of stability of the 

315 synchronized state may be valid for both smooth and non-smooth local dy- 

316 namics. Thus, results as those found in reference j23| are directly extended 

317 to piecewise-linear coupled maps. We have used general arguments valid for 

318 hyperbolic CML's so as to prove that the local transversal stability of the 

319 synchronized state actually implies the synchronization of all typical orbits. 



16| that the expression ([7]) is valid for 



320 Finally, since there is the conjecture 

321 non-hyperbolic systems, we conjecture that our results are valid for networks 

322 whose phase space has no structure hyperbolic. 
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Figure 1: (color online) Values of the coupling strength at the onset of transversal stability 
loss of the synchronized state, as a function of the effective coupling range. We used (3—1.1 
and (a) different lattice sizes; (b) different transient times Tw, for a fixed lattice size; (c) 
different distributions of initial conditions, for N = 257; (d) different initial conditions, for 
A'^ = 33 and a fixed transient time. The solid lines represent the analytical results from 
linear transversal stability of the synchronized state. 
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Figure 2: (color online) Line: decay in the number of trajectories in the unsynchronized 
state for N = 17 , ji = 1.10, a = 0.80 and e = 0.17 (slighty above the critical curve in 
Figure[T]). Histogram; synchronization times distribution for same lattice parameters. The 
exponential tail of such distribution is given by the escape rate of the saddle embedded in 
r. The average synchronization time, (n) — 98847, is aproximately given by jjj^f 
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Figure 3: (color online) (a) Series of a typical trajectory, x„, in the vicinity of a periodic 
orbit p„ G F. The trajectory follows the periodic orbit for a few periods, and then escapes 
from the saddle and reach the synchronized state. (b)The distance between the trajectory 
and the periodic orbit as a function of time n. Trajectory moves away from the UPO at 
a rate given by the unstable eigenvalue of the orbit. 
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Figure 4: (color online) (a) Density projection of the chaotic saddle that is contained in 
r, for typical trajectories. The color scale indicates the density: the blue color indicates 
regions most visited and the green color indicates the least visited, (b) Same as previous, 
but using all unstable orbits with the period less than or equal to 35. The simulations 
were performed with the following parameters: N — 17, f3 ^ 1.10, a = 0.80 and e = 0.17. 
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